Sonic Mach Cones Induced by Fast Partons in a Perturbative Quark-Gluon Plasma 
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We derive the space-time distribution of energy and momentum deposited by a fast parton travers- 
ing a weakly coupled quark-gluon plasma by treating the fast parton as the source of an external 
color field perturbing the medium. We then use our result as a source term for the linearized hy- 
drodynamical equations of the medium. We show that the solution contains a sonic Mach cone and 
a dissipative wake if the parton moves at a supersonic speed. 
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An interesting problem in the physics of deconfincd 
strongly interacting matter - the quark-gluon plasma - 
is to calculate the response of the medium to the passage 
of a fast parton, i. e., a parton with velocity approach- 
ing the speed of light. Fast partons are created experi- 
mentally in high-energy collisions of two nuclei when two 
energetic partons of the nuclear wave functions scatter 
at a large angle and acquire a large transverse momen- 
tum relative to the beam direction. If the interaction 
happens near the nuclear surface, one parton (commonly 
called the "trigger jet") rapidly leaves the medium and 
decays into a jet of hadrons which are cleanly observed in 
the detectors, while the second parton (commonly called 
the "back-jet") propagates with high velocity through 
the medium and deposits energy and momentum along 
the way in a process known as jet quenching (see, e. g., 
[H, HJ). Experimental measurements 0,0] of di-hadron 
correlation functions show a double peak structure in the 
back-jet distribution which is consistent with the forma- 
tion of a Mach cone-shaped emission pattern. 

The question of how the energy and momentum de- 
posited by the back-jet affects the bulk behavior of an 
evolving quark-gluon plasma has been examined by sev- 
eral authors (see, e. g., @, |, 0, I, I, M ED, d, E3). 
Using a schematic source in a two-dimensional hydrody- 
namics simulation of an expanding quark-gluon plasma, 
Chaudhuri and Heinz failed to find the formation of a 
Mach cone except for extreme values of the energy depo- 
sition [T3|. On the other hand, Casalderrey-Solana et al. 
Q showed that if one couples an appropriately chosen 
supersonic sound source to a linearized hydrodynamical 
equation one obtains a propagating Mach cone. However, 
the derivation of the source term expected from a super- 
sonic parton moving through a QCD plasma remained 
an open problem. 

Here we will derive this source term by calculating the 
effect of the color field of a fast parton on a perturbative 
quark-gluon plasma. We then solve the linearized hy- 
drodynamical equations for a static quark-gluon plasma 
and analyze how the perturbation propagates through 
the medium. For sound propagation we will use values of 
the sound attenuation length, r s , which lie in the range 
compatible with perturbation theory. These include the 



value obtained from the leading-order perturbative result 
for the shear viscosity rj [l4[ and the value for rj obtained 
by including three-body scattering processes in the trans- 
port cross section [l5|, [l6| . We note that the latter value 
is more consistent with shear viscosity allowed by the 
RHIC data [H • 

As explained above, we consider the fast parton as the 
source of an external color field that interacts with the 
medium. Asakawa et al. [H| showed that in the presence 
of soft color fields the color singlet parton distribution in 
a perturbative quark-gluon plasma obeys the following 
dissipative Vlasov-Boltzmann equation: 



^ + |-V x -V Pi Ai(p,*)V p , 



/(x,p,i) = C[/] 



(1) 

where /(x, p, t) is the ensemble averaged phase space dis- 
tribution of medium partons, p/E is the velocity of a 
parton with momentum p and energy E, and 



Dij(p,t)= f Fi{x.,t)Fj{x! ,t'), 

J — OO 



(2) 



where Fi (x, t) is the color Lorentz force on a medium 
particle: 



F,(x, t) = gQ a (t) (£7f(x, t) + (v x B)?(x, t)) 



(3) 



Here we consider the color fields in (|3|) to be generated 
by the fast parton. For a parton moving with velocity u 
with respect to the medium these fields, to lowest order in 
the coupling constant <?, have the Fourier representation 



E -(x,t) = ^| Uk e -i fc .*(""-kWa,-k-u) 
(27r) ci J k — lu 2 



B a (x,i) 



(2tt) 3 
(2^)3 



(4) 



d A k e 



. x (kx u)(S(w - k • u)) 



k 2 



(5) 

where fc M = (ui, k). In (|31[5]) we have neglected the screen- 
ing of the color field by the medium. We will incorporate 
the effect of color screening and short-distance quantum 
effects by appropriate infrared and ultraviolet cut-offs. 
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The hydrodynamical equations for the medium are ob- 
tained in the usual Chapman-Enskog approach by tak- 
ing moments of the evolution equation (fTJ). We assume 
that the medium is in local thermal equilibrium, with an 
energy-momentum tensor given by 



T" v = {e+p)u»u v ~pg^ - V {V^u w + V y u^ - d a u a ) 

(6) 

where = 7(1, v) denotes the medium four-velocity, 
= u^u" - g» u , V" = A» u d u , and T) is the shear 
viscosity. We evaluate each term by boosting to a frame 
co-moving with our volume element and then exploiting 
the assumption of local thermal equilibrium. Introducing 
the notation 

r = J ^(Vp.A^OV^p,*)), (7) 

we find that the resulting equations of motion for the 
medium evolution are 



r 



(8) 



where J u represents a source term due to the interaction 
of the medium with the passing fast parton. 

In order to solve the hydrodynamical equations ((5J), we 
assume that the energy and momentum density deposited 
by the parton is small compared to the equilibrium en- 
ergy density of the medium. This assumption permits us 
to linearize the hydrodynamical equations following the 
approach introduced in [6[. We write the perturbations 
of the energy-momentum tensor resulting from the source 
in an otherwise static medium as Se = ST oa and g with 
cji = ST ' 1 . The solutions of the linearized hydrodynami- 
cal equations can then be expressed in momentum space 
as 

ik,J L (k,uj) + J°(k,w)(iw- T s k 2 ) 



(5e(k,w) 

gi(k,w) 
gr(k,w) 



lu 2 — c 2 k 2 + iT s ojk 2 
iLvkJ L (k,u)) + ic 2 kJ°(k,w) 



CO 



2 — c 2 k 2 + iT s tok 2 



iJr(k, uj) 
■jj + \iT s k 2 



(9) 
(10) 
(11) 



where eo and po are the unperturbed energy density and 
pressure, respectively, c s denotes the speed of sound, and 
is the sound attenuation length. 
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" s — 3(e +Po) 3s T 

In the above equations we have divided the source and 
perturbed momentum density vectors into transverse and 
longitudinal parts: g = gL + gr and J = k Jl + Jt, with 
k denoting the unit vector in the direction of k. We note 
that (fTTj) is the diffusion equation and does not describe 
sound propagation. If a Mach cone appears it will be 
found in the dynamics of eqs. (f9l fT0|) . 

In order to proceed further we need to evaluate the 
source term ([7])- We consider a thermal plasma of mass- 
less gluons with the unperturbed distribution 



/o(x,p,i) 



2{N 2 - 1) 
ePv° - 1 



(12) 



where p° = E = |p| and N c = 3 is the number of col- 
ors. Ignoring the time dependence of Q a (t) in ([3]) we 
eventually find that 



J°(x,t) 



(2tt) s 



g- / d 4 fcd 4 /c'e l ( k+k '> x - l ^+"') t ; 



rfv 



(v-E a (fc'))(v-E a (fc)) 
47r(w' - k' • v + ie) 



(13) 



and 



J fc (x,£) = 



(2tt) s 



d 4 fcd 4 fcV (k+k,) - x - 4(w+cy) *x 



r/v 



(v-E°(fcQ)(^(fc) + (vxB)g(fc)) 
4it(uj' — k' • v + ie) 



(14) 



where the fields E a (/c) and B a (k) are given by ([4]) and 
((5|), respectively. The Debye mass for gluons is given 
by mo = gT. If one includes massless quarks in the 
medium the expressions (|13l and I14[) remain unchanged 
except that the quark contribution to the Debye mass 
needs to be included in tod ■ The longitudinal part of J 
is obtained by multiplication with k: Jl — k • J. The 
source distributions can be analytically evaluated in the 
ultra-relativistic limit 7 = (1 — it 2 ) -1 / 2 3> 1. After a 
lengthy calculation one obtains for a parton moving in 
the positive z-direction with velocity u: 



~fu{z — ut) 



I , (15) 

z — ut 



J°(p,Z,t) = d(p,z,t)ju 2 ( 

J z (p,z,t) = uJ°(p, z,t) — d(p, z,t)u 2 , (16) 

P 



3j_(p,z,t) 



—d(p, z, t)u 2 —^, 
P 



(17) 



where we introduced the notation xj^ = (x,y) and p 
|xj_|. The function d is given by 



d(p,z,t) 



OisQlm 2 ^ 



8tt[p 2 + 7 2 (z-irf) 2 ] 3 / 2 ' 



(18) 



encoding the Lorentz contracted field configuration of the 
fast parton. A detailed derivation of the source term, 
including color screening by the medium, is presented in 
[20I ]. Integrating over all space and introducing infrared 
and ultraviolet cut-offs for the p-integration, we obtain 
the differential energy loss 



Pmin 



(19) 



For p max = l/m D and p min = l/(2y/E p T), where E p is 
the energy of the fast parton, one recovers the standard 
leading-logarithmic result [lj| for elastic energy loss in 
the medium. A stronger than logarithmic dependence of 
the energy loss on E p , found in some other approaches, 
may be effectively described by a parametric dependence 
of a s Q 2 on E p . 
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Since hydrodynamics is only valid at distances large 
compared to the mean free path, and in a weakly cou- 
pled plasma the mean free path is parametrically large 
compared to the color screening length, the source term 
generated by an energetic parton is, in first approxima- 
tion, point-like. In this spirit, the source term (jl5til8[) 
derived here can be thought of as a sophisticated repre- 
sentation of a delta function. However, it also contains, 
in the sense of effective field theory, an infinite series 
of higher derivative contributions, which produce sub- 
leading corrections. 

At leading order in the gradient expansion, both the 
amplitude and structure of the Mach cone produced by 
the projectile are uniquely determined by the rate of elas- 



tic energy loss together with the value of viscosity. The 
source term derived and used here provides a model of the 
sub-leading momentum dependence of the source term 
for the special case of a weakly coupled plasma, which 
approximately captures the physics at the Debye screen- 
ing scale, but ignores physics at longer sub-hydrodynamic 
length scales |21| . 

To calculate the propagating disturbance of the 
medium by the projectile we insert (|13p and (|14jl into 
eqs. (|51-|lip. For a relativistic parton, propagating with 
velocity u = (0, 0, u) and position x = (0, 0, ut) we ob- 
tain the following expressions for the quantities (j9l- [TTj) . 
Fourier transformed back to space-time: 



_ mVQ|m| r ik . (x _ xVittfezt / -k-^ + (2uk z+i T s k 2 )(p' 1 -z' 1 2 u) \ 

1 ' ' ~ 327r 2 (27r) 3 (u 2 - c 2 ) J aK ° X [ (k 2 - \ 2 k 2 + ia)p'(p' 2 + 7 V 2 )3/ 2 J ' {M > 



327r 2 (27r)3( u 2 _ c 2) J I fe2 ( fc 2 _ A 2 fc 2 + ia y{ p ii + ^2^3/2 



, * 1 [ M , , ik.(x-x')-i U fc t J(x')fc 2 -k(k- J(x')) 



(22) 



where A 2 = c 2 J{u 2 - c 2 ), a = T s u(\ 2 /c 2 )k z {k$ + A: 2 ), 
and Q 2 is the Casimir operator for the color charge 
of the projectile, which is 4/3 for a quark and 3 for 
a gluon. In our evaluation of (|20j|22p we perform five 
of the six integrals analytically leaving the final inte- 
gral over kr = {k 2 + k 2 ) 1 / 2 to be done numerically. 
In (|2"D1 |2"T} we perform the integral over k z using con- 
tour integration and find poles at (among other places) 
k z = ±(k 2 ,\ 2 =p ilcl) 1 / 2 , where \a\ is itself a function of 
k z . We approximate a(k z ) by o^ifcrA) when evaluating 
the residues at these poles, which is permitted at mo- 
mentum scales for which the sound attenuation is small 
(kr <C c 2 /r s ). We note that our expression contains the 
term (p' 2 +7 2 z' 2 ) -3 / 2 which is (up to appropriate normal- 
ization) a nascent delta function in z' in the limit 7 — > 00. 
With this insight, and noting that we are working in the 
ultrarelativistic limit, we expand the z 1 dependence of the 
exponent to first order to obtain analytically manageable 
expressions. An ultraviolet cut-off of the order of T s in 
the p' integration and an appropriate infrared cut-off in 
the z' integration are used to regularize the logarithmic 
divergences when necessary. 

We now present and discuss the numerical results ob- 
tained from our evaluation of (|20H22p for the case of 
a gluon moving along the positive z axis with velocity 
u = 0.99955c (7 ~ 33). As mentioned in the beginning, 



we will use two different values of the sound attenua- 
tion length, r s , which have been calculated perturba- 
tively and compare the results. The first value is de- 
termined from the leading order result for rj/s obtained 
by Arnold et al. [Til ]. This result includes only binary 
(2^2 and 1 — > 2) processes and gives 77/s = 0.48 for 
a gluonic plasma with a s = 0.3. The second value for 
77/s is the one obtained by Xu and Greiner by going be- 
yond leading order in the diluteness of the medium [22| . 
Their calculation, including 3 <-> 2 scattering processes, 
found 77/s = 0.13 for a gluonic plasma with a s = 0.3 
[l5| . Even smaller values of 77/s may be compatible with 
the assumption of a perturbative medium, if the viscos- 
ity is lowered by anomalous contributions [23| . Finally, 
we have chosen T = 350 MeV as the temperature of the 
medium and use c s = c/V3 for the speed of sound. 

Figures Q] and [2] show the scaled energy density (left) 
and momentum density (right) deposited in the medium 
by the projectile along the plane (x, 0, z — ut) for rj/s 
0.1 .') and 77/s = 0.48, respectively. In both cases the 
energy density wave excited in the medium is seen to de- 
velop the shape of a Mach cone whose intensity is peaked 
near the source. At growing distance from the source the 
Mach cone slowly weakens and broadens. As one would 
expect, the cone structure is more pronounced for the 
smaller value of 77/s. It is instructive to compare the 



4 




FIG. 1: (Color online) Plots of the scaled perturbed energy density |x|<5e(x)/(mpT) (left), and momentum density 
|x|g(x) /(m^T) (right), excited by a gluon moving along the positive z axis at position ut and speed u = 0.99955c for i]/s = 0.13. 
The results, which are cylindrically symmetric around the z axis, are shown on the plane x = (a;,0, z — ut). The energy and 
momentum densities have been multiplied by the distance |x| from the source to compensate for the geometric dilution effect 
of the cone. The values chosen for the parameters mo, ot a , T, c 3 , and F s are disussed in the text. 




FIG. 2: (Color online) The same as in Figure [T] but with r\js — 0.48. Just as in Figure [T] the energy density wave excited in 
the medium develops the shape of a Mach cone. Here, however, the distribution is more damped and spread out. 



wavelength below which sound propagation is strongly 
damped, A c = 2t:T s /c s , with the extent of the sound 
source, m^ 1 sa 0.3 fm. For rj/s = 0.48 we have A c /4 m 1.0 
fm, while for rj/s = 0.13 we find A c /4 « 0.25 fm, suggest- 
ing that the coupling of the source term to the sound 
mode becomes increasingly ineffective for rj/s 0.15. 

We next consider the velocity fields induced by the 
moving source which are given by v = g/eo where 



eq 10 GeV/fm 3 is the energy density of a plasma of 
massless gluons at T = 350 MeV. We plot the induced 
velocity flow in Fig. [3] for both the contribution from the 
sound equation ()21|) and that from the diffusion equa- 
tion |22|) for the case of 77/s = 0.13. The velocity flow 
obtained from the sound equation has a similar struc- 
ture as the perturbed energy density with a well defined 
Mach cone trailing behind the moving parton. As men- 
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tioncd before, the velocity flow obtained from the diffu- 
sion equation does not describe sound propagation but 
instead describes the source pulling matter along in its 
wake. For our choice of parameters the collective flow 
induced by the diffusive wake has a longitudinal velocity 
of approximately O.f 5c at a distance of 4 fm behind the 
source and a much smaller transverse velocity. 

In summary, we have derived the pattern of energy 
and momentum deposition by a fast parton traversing a 
quark-gluon plasma in perturbative QCD. Our result de- 
pends on the following parameters: the source strength 



Q 2 p , the Debye mass mo, the sound velocity c s , and 
the sound attenuation length, T s . Treating the propaga- 
tion of the disturbance created by the projectile in lin- 
earized viscous hydrodynamics we have shown that the 
fast moving parton excites a sonic Mach cone and a dif- 
fusive wake. The intensity of the Mach cone was found 
to decrease with growing kinematic viscosity. 
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An interesting problem in the physics of deconfincd 
strongly interacting matter - the quark-gluon plasma - 
is to calculate the response of the medium to the passage 
of a fast parton, i. e., a parton with velocity approach- 
ing the speed of light. Fast partons are created experi- 
mentally in high-energy collisions of two nuclei when two 
energetic partons of the nuclear wave functions scatter 
at a large angle and acquire a large transverse momen- 
tum relative to the beam direction. If the interaction 
happens near the nuclear surface, one parton (commonly 
called the "trigger jet") rapidly leaves the medium and 
decays into a jet of hadrons which are cleanly observed in 
the detectors, while the second parton (commonly called 
the "back-jet") propagates with high velocity through 
the medium and deposits energy and momentum along 
the way in a process known as jet quenching (see, e. g., 
[H, 0)- Experimental measurements 0,0] of di-hadron 
correlation functions show a double peak structure in the 
back-jet distribution which is consistent with the forma- 
tion of a Mach cone-shaped emission pattern. 

The question of how the energy and momentum de- 
posited by the back-jet affects the bulk behavior of an 
evolving quark-gluon plasma has been examined by sev- 
eral authors (see, e. g., @, |, 0, I, I, M ED, d, E3). 
Using a schematic source in a two-dimensional hydrody- 
namics simulation of an expanding quark-gluon plasma, 
Chaudhuri and Heinz failed to find the formation of a 
Mach cone except for extreme values of the energy depo- 
sition [T3|. On the other hand, Casalderrey-Solana et al. 
Q showed that if one couples an appropriately chosen 
supersonic sound source to a linearized hydrodynamical 
evolution one obtains a propagating Mach cone. In the 
above cases the explicit form of the source term expected 
from a supersonic parton moving through a QCD plasma 
was not addressed. 

Here we will derive the magnitude and shape of this 
source term, valid at distances below the Debye screen- 
ing length, by calculating the effect of the color field of 
a fast parton on a perturbative quark-gluon plasma. We 
then solve the linearized hydrodynamical equations for 
a static quark-gluon plasma and analyze how the per- 
turbation propagates through the medium. For sound 
propagation we will use values of the sound attenuation 



length, r s , which lie in the range compatible with pertur- 
bation theory. These include the value obtained from the 
leading-order perturbative result for the shear viscosity 
77 [3] and the value for r\ obtained by including three- 
bodyscattering processes in the transport cross section 
IH 16| . We note that the latter value is more consistent 



17] 



with shear viscosity allowed by the RHIC data 

As explained above, we consider the fast parton as the 
source of an external color field that interacts with the 
medium. Asakawa et al. fl8l | showed that in the presence 
of soft color fields the color singlet parton distribution in 
a perturbative quark-gluon plasma obeys the following 
dissipative Vlasov-Boltzmann equation: 



d 

^ + |-V a ,-V ft D«(p,t)V ft 



/(x,p,i) = C[/] 



(1) 

where /(x, p, t) is the ensemble averaged phase space dis- 
tribution of medium partons, p/E is the velocity of a 
parton with momentum p and energy E, and 



D ij (p,t)= j dt'Fi^F^J), 

J —OO 



(2) 



where Fi (x, t) is the color Lorentz force on a medium 
particle: 

Fi(x,t) = gQ a (t) (E?( x ,t) + (v x B)?(x,i)) . (3) 

Here we consider the color fields in to be generated 
by the fast parton. For a parton moving with velocity u 
with respect to the medium these fields, to lowest order in 
the coupling constant g, have the Fourier representation 



(2tt) 3 



k 2 



(4) 



{2-K) i J k l — u z 

(5) 

where fc M = (ui, k). In (jH[5|) we have neglected the screen- 
ing of the color field by the medium. We will incorporate 
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the effect of color screening and short-distance quantum 
effects by appropriate infrared and ultraviolet cut-offs. 

The hydrodynamical equations for the medium are ob- 
tained in the usual Chapman-Enskog approach by tak- 
ing moments of the evolution equation ([I]). We assume 
that the medium is in local thermal equilibrium, with an 
energy-momentum tensor given by 

T» v = {e+p)u^u v -pg^ -r}(Wu" + W--A^d a u Q ) 

3 

(6) 

where — 7(1, v) denotes the medium four- velocity, 
W = vTu" - g^, V = A» v d u , and r) is the shear 
viscosity. We evaluate each term by boosting to a frame 
co-moving with our volume element and then exploiting 
the assumption of local thermal equilibrium. Introducing 
the notation 



r = 



dp p h 
~{2nf 



(V Pi Ai(p,*)V P3 ./(x,p,t)) J (7) 



we find that the resulting equations of motion for the 
medium evolution are 



r 



(8) 



where J u represents a source term due to the interaction 
of the medium with the passing fast parton. 

In order to solve the hydrodynamical equations |8[). we 
assume that the energy and momentum density deposited 
by the parton is small compared to the equilibrium en- 
ergy density of the medium. This assumption permits us 
to linearize the hydrodynamical equations following the 
approach introduced in [||. We write the perturbations 
of the energy-momentum tensor resulting from the source 
in an otherwise static medium as 8e = ST m and g with 
gi = ST ' 1 . The solutions of the linearized hydrodynami- 
cal equations can then be expressed in momentum space 
as 



<5e(k, uj) 
gi(k,w) 
gr(k,w) 



ik J L (k, uj) + J°(k, uj)(iuj - T s k 2 ) 

uj 2 - c 2 k 2 + iT s ujk 2 
iujkJ L (k,u)) + ic 2 s kJ°(k,uj) 
uj 2 — c 2 k 2 + iV s ujk 2 
iJy(k, U)) 



■iT,k 2 



. (9) 
(10) 
(11) 



where eo and po are the unperturbed energy density and 
pressure, respectively, c s denotes the speed of sound, and 
r s = 57 4 '' — r = -^Kp is the sound attenuation length. 

* 3(eo+Po) 3s T ° 

In the above equations we have divided the source and 
perturbed momentum density vectors into transverse and 
longitudinal parts: g = gL + gT and J = k Jl + 3t, with 
k denoting the unit vector in the direction of k. We note 
that (fTTj) is the diffusion equation and does not describe 
sound propagation. If a Mach cone appears it will be 
found in the dynamics of eqs. (f^l fTOj) . 

In order to proceed further we need to evaluate the 
source term ([7])- We consider a thermal plasma of mass- 



less gluons with the unperturbed distribution 

2(N 2 C 1) 



/o(x,p,i) 



2 /3p° _ 1 



(12) 



where p° = E = |p| and N c = 3 is the number of col- 
ors. Ignoring the time dependence of Q a (t) in ([3]) we 
eventually find that 



J°(x,*) = 



(2tt) s 



d 4 fcd 4 /fcV (k+k,) - x - l(wW)t 



(v-E a (fc'))(v-E a (fc)) 
Att(uj' — k' ■ v + it) 



(13) 



and 



J fc (x,i) = 



(2tt) 6 



d 4 fcd 4 fcV( k+k ')- x -^+^'x 



r/v 



(v-E°(fc'))(£°(fc) + (vxB)°(fc)) 
47r(w' — k'-v + ie) 



(14) 



where the fields ~E a (k) and B a (/e) are given by ([?]) and 
([5]), respectively. The Debye mass for gluons is given 
by TUj) = gT. If one includes massless quarks in the 
medium the expressions ((TBI and PH]) remain unchanged 
except that the quark contribution to the Debye mass 
needs to be included in mo- The longitudinal part of J 
is obtained by multiplication with k: Jl = k ■ J. The 
source distributions can be analytically evaluated in the 
ultra-relativistic limit 7 = (1 — it 2 ) -1 / 2 3> 1. After a 
lengthy calculation one obtains for a parton moving in 
the positive z-direction with velocity u: 

J°(p,z,t) = d{p,z,t) lu 2 fl- T^r^ V (15) 



, (16) 



J z (p,z,t) = uJ (p, z,t) — d(p, z,t)u 



. z — ut 



3±(p,z,t) 



-d(p, z, t)u 



2 X -L 



(17) 



where we introduced the notation xj^ = (x,y) and p 
|xj_|. The function d is given by 



d(p,z,t) 



asQlm 2 ^ 



8n[p 2 + 1 2 {z-ut) 2 } 3 / 2 ' 



(18) 



encoding the Lorentz contracted field configuration of the 
fast parton. A detailed derivation of the source term, 
including color screening by the medium, is presented in 
[20j . Integrating over all space and introducing infrared 
and ultraviolet cut-offs for the p-integration, we obtain 
the differential energy loss 



f ,3 T o / \ a s ^2 2 l Pmax 

— = / d 6 xr{x) = —Q^mfyln . 

dX I 2 y Pmin 



(19) 



For p max = l/m D and p min = l/(2y/E p T), where E p is 
the energy of the fast parton, one recovers the standard 
leading-logarithmic result [lj| for elastic energy loss in 



3 



the medium. A stronger than logarithmic dependence of 
the energy loss on E p , found in some other approaches, 
may be effectively described by a parametric dependence 
of a s Q 2 on E p . 

Since hydrodynamics is only valid at distances large 
compared to the mean free path, and in a weakly cou- 
pled plasma the mean free path is parametrically large 
compared to the color screening length, the source term 
generated by an energetic parton is, in first approxima- 
tion, point-like. In this spirit, the source term (|15H18[) 
derived here can be thought of as a sophisticated repre- 
sentation of a delta function. However, it also contains, 
in the sense of effective field theory, an infinite series 
of higher derivative contributions, which produce sub- 
leading corrections. 

At leading order in the gradient expansion, both the 



amplitude and structure of the Mach cone produced by 
the projectile arc uniquely determined by the rate of clas- 
tic energy loss together with the value of viscosity. The 
source term derived and used here provides a model of the 
sub-leading momentum dependence of the source term 
for the special case of a weakly coupled plasma, which 
approximately captures the physics at the Debye screen- 
ing scale, but ignores physics at longer sub-hydrodynamic 
length scales |2l| . 

To calculate the propagating disturbance of the 
medium by the projectile we insert (TT3"|) and (TT4")) into 
eqs. (|9l-|lip. For a relativistic parton, propagating with 
velocity u = (0, 0, u) and position x = (0, 0, ut) we ob- 
tain the following expressions for the quantities (Pl- fTTj) . 
Fourier transformed back to space-time: 



iu 2 g 2 Q 2 p m 2 D f , lk . (x _ x <)_ mfc t /-k ■ x' + (2uk z + iT s k 2 ) (p'j - z'^ 2 u)\ 
= _______________ y rfkrfx e - xx «-l ___________________ I (20) 



iu 2 g 2 Q 2 p m 2 D r ik . (x _ x0 _ iMiM ( ~uk z k ■ x' + ((uk z ) 2 + c 2 k 2 ) {pfj " \ , 2n 

32^ 2 (2^) 3 (u 2 -c 2 ) J I fc 2 (fc 2 - A 2 fc 2 + ia)p'{p' 2 + 1 2 z' 2 ) 3 / 2 J 1 ' 

gT(M) . ^/^M.-.,-.,, ^^) (22) 



where A 2 = c 2 /(u 2 - c 2 ), a = Y s u{\ 2 j c 2 s )k z (k\ + k 2 ), 
and Q 2 is the Casimir operator for the color charge 
of the projectile, which is 4/3 for a quark and 3 for 
a gluon. In our evaluation of (|20j|22p we perform five 
of the six integrals analytically leaving the final inte- 
gral over kr = {k 2 + k 2 ) 1 / 2 to be done numerically. 
In (|2"D1 l_Tj) we perform the integral over k z using con- 
tour integration and find poles at (among other places) 
k z = ±(fc^A 2 =F ilcl) 1 / 2 , where \a\ is itself a function of 
k z . We approximate a(k z ) by a^ifcrA) when evaluating 
the residues at these poles, which is permitted at mo- 
mentum scales for which the sound attenuation is small 
(kr <C c 2 /r s ). We note that our expression contains the 
term (p' 2 +7 2 z' 2 ) -3 / 2 which is (up to appropriate normal- 
ization) a nascent delta function in z' in the limit 7 — ► 00. 
With this insight, and noting that we are working in the 
ultrarelativistic limit, we expand the z' dependence of the 
exponent to first order to obtain analytically manageable 
expressions. An ultraviolet cut-off of the order of T s in 
the p' integration and an appropriate infrared cut-off in 
the z' integration are used to regularize the logarithmic 
divergences when necessary. 

We now present and discuss the numerical results ob- 
tained from our evaluation of (|20H22p for the case of 



a gluon moving along the positive z axis with velocity 
u = 0.99955c (7 w 33). As mentioned in the beginning, 
we will use two different values of the sound attenua- 
tion length, T s , which have been calculated pcrturba- 
tively and compare the results. The first value is de- 
termined from the leading order result for r//s obtained 
by Arnold et al. [TJ]. This result includes only binary 
(2^2 and 1 — * 2) processes and gives r\js = 0.48 for 
a gluonic plasma with a s = 0.3. The second value for 
r\j s is the one obtained by Xu and Greiner by going be- 
yond leading order in the diluteness of the medium [_| . 
Their calculation, including 3^2 scattering processes, 
found r\j s = 0.13 for a gluonic plasma with a s = 0.3 
fl5| . Even smaller values of ?]/s may be compatible with 
the assumption of a perturbative medium, if the viscos- 
ity is lowered by anomalous contributions [__| . Finally, 
we have chosen T = 350 MeV as the temperature of the 
medium and use c s = c/V3 for the speed of sound. 

Figures Q] and [2] show the scaled energy density (left) 
and momentum density (right) deposited in the medium 
by the projectile along the plane (x, 0, z — ut) for r\j s = 
0.13 and r\js = 0.48, respectively. In both cases the 
energy density wave excited in the medium is seen to de- 
velop the shape of a Mach cone whose intensity is peaked 



4 




FIG. 1: (Color online) Plots of the scaled perturbed energy density |x|<5e(x)/(mpT) (left), and momentum density 
|x|g(x) /(m^T) (right), excited by a gluon moving along the positive z axis at position ut and speed u = 0.99955c for i]/s = 0.13. 
The results, which are cylindrically symmetric around the z axis, are shown on the plane x = (a;,0, z — ut). The energy and 
momentum densities have been multiplied by the distance |x| from the source to compensate for the geometric dilution effect 
of the cone. The values chosen for the parameters mo, ot a , T, c 3 , and F s are disussed in the text. 




FIG. 2: (Color online) The same as in Figure [T] but with r\js — 0.48. Just as in Figure [T] the energy density wave excited in 
the medium develops the shape of a Mach cone. Here, however, the distribution is more damped and spread out. 



near the source. At growing distance from the source the 
Mach cone slowly weakens and broadens. As one would 
expect, the cone structure is more pronounced for the 
smaller value of rj/s. It is instructive to compare the 
wavelength below which sound propagation is strongly 
damped, A c = 2ixT s /c s , with the extent of the sound 
source, m^ 1 ps 0.3 fm. For 77/s = 0.48 we have A c /4 w 1.0 
fm, while for rj/s = 0.13 we find A c /4 f» 0.25 fm, suggest- 



ing that the coupling of the source term to the sound 
mode becomes increasingly ineffective for rj/s S> 0.15. 

We next consider the velocity fields induced by the 
moving source which are given by v = g/eo where 
eo ~ 10 GeV/fm 3 is the energy density of a plasma of 
massless gluons at T = 350 MeV. We plot the induced 
velocity flow in Fig. [3] for both the contribution from the 
sound equation (|21[) and that from the diffusion equa- 
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tion (|22|) for the case of r\js — 0.13. The velocity flow 
obtained from the sound equation has a similar struc- 
ture as the perturbed energy density with a well defined 
Mach cone trailing behind the moving parton. As men- 
tioned before, the velocity flow obtained from the diffu- 
sion equation does not describe sound propagation but 
instead describes the source pulling matter along in its 
wake. For our choice of parameters the collective flow 
induced by the diffusive wake has a longitudinal velocity 
of approximately 0.15c at a distance of 4 fm behind the 
source and a much smaller transverse velocity. 

In summary, we have derived the pattern of energy 
and momentum deposition by a fast parton traversing a 



quark-gluon plasma in pcrturbative QCD. Our result de- 
pends on the following parameters: the source strength 
a s Qp, the Debye mass mo, the sound velocity c s , and 
the sound attenuation length, T s . Treating the propaga- 
tion of the disturbance created by the projectile in lin- 
earized viscous hydrodynamics we have shown that the 
fast moving parton excites a sonic Mach cone and a dif- 
fusive wake. The intensity of the Mach cone was found 
to decrease with growing kinematic viscosity. 
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FIG. 3: (Color online) Plots of the velocity field for the same scenario as in Figure [T] The two lower plots indicate the parallel 
(i.e. moving with the source) components of the induced velocity flows for both the sound equation (|2ip and the diffusion 
equation (|22[) while the two upper plots are for the perpendicular components. The axes are plotted in units of inverse GeV. 



